Impact of Sequential Disorder on Scaling Behavior of Boarding Time 
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Airplane boarding process is an example where disorder properties of the system are relevant to 
emergence of universality classes. Based on a simple model, we present a systematic analysis of finite- 
size effects in boarding time, and propose a comprehensive view of the role of sequential disorder 
in the scaling behavior of boarding time against the plane size. Using numerical simulations and 
mathematical arguments, we find how the scaling behavior depends on the number of seat columns 
and the range of sequential disorder. Our results show that new scaling exponents can arise as 
disorder is localized to varying extents. 

PACS numbers: 89.75.Da, 89.40.Dd, 02.50.-r, 05.40.-a 
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How disorder properties affect the scaling behavior of 
a characteristic time scale is a problem studied in vari- 
ous areas of physics. For example, dynamic exponents of 
surface growth models [ij are known to depend on the 
geometry of substrate and the type of disorder. In com- 
plex networks, consensus time of opinion dynamics @] 
and first-passage time of random walk Q scale differently 
with the system size depending on the level of structural 
heterogeneity. Those studies reveal the subtle interplay 
between dynamics and disorder that gives rise to differ- 
ent universality classes of scaling behaviors, each of which 
is succinctly characterized by its scaling exponents. Ex- 
panding the list of those universality classes and clarify- 
ing their origins has been established as a general frame- 
work frequenty employed by physicists. 

Airplane boarding process provides another interesting 
example where this framework can be applied. The av- 
erage time required for all passengers to get seated, or 
average boarding time (T) , may scale differently with the 
plane size N depending on sequential disorder of passen- 
gers, which is controlled by the airline's boarding policy 
and each individual's gate arrival time. However, most 
studies of boarding process were limited to the practi- 
cal problem of reducing boarding time at a fixed plane 
size (^-fioj. which provides only fragmented knowledge 
about performances of different boarding policies in pre- 
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FIG. 1: Snapshots of boarding in an airplane with TV = 3 rows 
and c — 2 columns of seats, where each passenger (circle) 
is labeled by the row index of its seat. At each time step, 
passengers are waiting in line (white), sitting down (black), 
or staying seated (gray). 



scribed situations. To understand the general nature of 
boarding process, we should pay attention to its scal- 
ing properties, which came to be studied only recently 
with inconsistent conclusions. While analytical studies 
by Bachmat et al. [Ill4l4| supported (r)^ iV^/^, a nu- 
merical study by Frette and Hemmer [l5| based on a 
simple model (see Fig. [1]) reported {T)^ ]\f0.69^ r^j^^^ 
was again contradicted by Bernstein [16|, who observed 
that the same model yields (T) N^^^ after filtering out 
finite-size effect. 

In this Letter, we resolve the conflict between differ- 
ent exponents by showing that T ^ N^^"^ is the correct 
description of the previous results if we carefully inter- 

"l. We 
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pret the relevant mathematical theorems 
also propose that different scaling behaviors are possi- 
ble if the range of disorder in the passenger sequence is 
localized to varying extents, which gives a more compre- 
hensive picture of the problem. To support these claims, 
we generalize the boarding model proposed by [l^, in- 
corporating additional parameters for the structure of the 
airplane and disorder in the passenger sequence. 

Model — We consider an airplane with N rows and 
c columns of seats along a one-dimensional aisle. Each 
seat is labeled with an index r if it is in the r-th row 
from the front. The aisle is discretized into N sites, each 
of which cannot hold more than one passenger at once. 
At t = 0, cN passengers enter the plane from the front 
in a sequence {ri, r2, . . . , r^jvlj where denotes the seat 
index of the i-th passenger to enter. Passengers' posi- 
tions arc synchronously updated. Every passenger moves 
along the aisle front-to-back, instantly crossing successive 
empty sites until blocked by an occupied site or reaching 
the row of its assigned seat. If the latter is the case, sit- 
ting down takes one time step for every passenger. Con- 
tinuing the process, all passengers would get seated at 
t = T, which we call boarding time (see Fig. [T]). 

Since the boarding process is deterministic in the 
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FIG. 2: Examples of the passenger sequence {r^}. As N ^ oo, 
(a) fixed N/m means globalized disorder; (b) fixed m means 
localized disorder; (c) Gaussian sliufHing gives either purely 
globalized or purely localized disorder; (d) uniform shuffling 
allows coexistence of both kinds of disorder. We set A'' = f 0* 
in (c) and (d). 



model, T is determined by the initial passenger sequence. 
Lack of strict control over the sequence implies that it 
is more meaningful to deal with (T) for an ensemble of 
sequences rather than T for a particular sequence [l9| . 
As one can understand from Fig. [21 such an ensem- 
ble can be represented as a probability density function 
(PDF) p{i/N,r/N), where i/N and r/N denote sequen- 
tial and seat indices of passengers normalized by N, re- 
spectively . This representation enables us to uti- 
lize the following mathematical theorem [TJ [S] that has 
profound implications for the scaling behavior of (T). 

Theorem — If (xq, ya), a = I, . . . , N are pairs of real 
numbers with < Xa < 1 and < < 1, we say that 
a subsequence {(x^^ , j/ij, . . . , (x^, , j/i,)} is an increasing 
subsequence if 



< Xi-^^ and yi. < y^.^, for j 
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where ij is a sequence of non-repeated indices between 1 
and A''. If the pairs {xa,ya) are generated from a finite 
PDF p{x,y), the length of the longest increasing subse- 
quence asymptotically scales as A^/^. 

The increasing subsequence in the theorem can be 
translated as the blocking subsequence of passengers, in 
which one passenger blocks the next one if the latter can- 
not reach its seat unless the former is seated H, IT^ . 



Since the length of the longest blocking subsequence is 
equal to boarding time T, the theorem indicates that T 



(and (T) as well) scales as N^/"^ if the PDF remains finite 
in the asymptotic limit N oo. The finitude of the PDF 
is equivalent to the broad range of disorder in the pas- 
senger sequence, which is comparable with the plane size 
and might be described as "globalized" . In other words, 
(T) ^ N^/"^ is guaranteed if sequential disorder is global- 
ized. Meanwhile, different scaling behaviors might arise 
if some part of the PDF diverges due to the presence of 
"localized" disorder in the passenger sequence. 

The range of sequential disorder is primarily deter- 
mined by boarding policies. A boarding policy is a set 
of constraints imposed on the initial passenger sequence. 
Here we focus on back-to-front policies with equal-sized 
boarding groups. According to such policies, passen- 
gers are divided into boarding groups of m rows each 
(m is a divisor of A^), so that the n-th group to enter 
the plane consists of passengers whose seat indices sat- 
isfy A — nm + 1 < < A — (ri — l)m. Consequently, 
the passenger sequence is sorted back-to-front on a scale 
larger than groups, but remains randomized within each 
group. The range of sequential disorder is given by to. 
If N/m (m) is fixed as A^ — cxd, sequential disorder is 
purely globalized (localized), as illustrated in Fig. [UJa) 
and (b), respectively. 

Meanwhile, some passengers would naturally deviate 
from the boarding policy due to their impunctuality. 
This justifies the incorporation of arrival time fluctua- 
tions as another determinant of sequential disorder. We 
consider Gaussian shuffling as a model of such fluctua- 
tions, which is defined as follows: we add an i.i.d. Gaus- 
sian random number rji of zero mean and of variance to 
the sequential index i of each passenger, and then sort the 
sequence in the increasing order of i + rji (see Fig. [HJc)). 
As long as to is finite, Gaussian shuffling makes sequen- 
tial disorder purely localized (globalized) if a {a/N) is 
fixed in the asymptotic limit. Alternatively, we can also 
consider uniform shuffling, in which every passenger is 
randomly relocated in the sequence with probability p 
(see Fig. [l{d)). If to is finite, uniform shuffling allows 
localized and globalized disorder to coexist until the se- 
quence is completely randomized at p = 1. 

Now that we have a collection of parameters control- 
ling the range of sequential disorder, we are ready to test 
numerically whether (T) N^^^ is true for purely glob- 
alized disorder, and to examine unkown scaling behaviors 
in the presence of localized disorder. We use the method 
of effective scaling exponent 0| to clarify the asymp- 
totic scaling, which is effectively the same as the method 
adopted by [16=|. 

Single-column airplane — We first consider a fully re- 
served plane with only a single column (c = 1) of seats. 
Although a single-column plane hardly exists in reality, 
it is not merely a theoretical construct. For example, 
consider a boarding policy that fills the window seats 
first and the aisle seats later. The policy separates the 
passengers to be seated in the same row very far from 
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FIG. 3: (Color online) Asymptotic 
scaling behaviors of (T) in a fully 
reserved plane. When c = 1, from 
bottom to top, the curves represent 
(a) m = 2, 4, 8, N/A, N/2, N with- 
out shuffling, and m — 10 with (b) 
a/N = 0.001,0.01,0.1,1, (c) a = 
1, 10, 100 (inset: semi-log plots of 
the main), (d) p = 0,0.01,0.1,1. 
When c — 2, from top to bottom, 
(e) m = 1,4, 16,iV/16,Af/4,iV 
without shuffling, and m = 5 
with (f) a/N = 0.001,0.01,0.1, 
(g) a = 1, 10, 100, 1000, (h) p = 
0, 0.5, 0.9, 1. Each data point is av- 
eraged over at least 100 samples. 



each other, so there is effectively no interference between 
them. In such a case, even a multi-column plane can be 
regarded as a single-column plane. 

The scaling behaviors of (T) for c = 1 are shown in 
Fig. [2Ia)~(d). Without arrival time fluctuations, (T) 
scales as TV^/^ if N/ m is fixed, while it saturates to m if 
m is fixed (see Fig. EJa)). Gaussian shuffling at finite m 
produces similar but slightly different scalings: the N^^'^ 
scaling for fixed a/N (see Fig.[3jb)) and the In TV scaling 
for fixed a (see Fig. [Sjc)). Since both saturation and log 
scaling are slower than any algebraic scaling, they can 
be collectively labeled as N'^ scalings. Thus, our results 
confirm (T) ^ N^^^ for purely globalized disorder, while 
finding (T) ^ N'^ for purely localized disorder. When 
those two kinds of disorder coexist, the effect of localized 
disorder predominates, as implied by the A^-^/^ scaling 
observed whenever p > (see Fig. IH^d)). 

We provide an argument for the origin of this new scal- 
ing in the special case when there are no arrival time 
fluctuations. As iV oo, the passenger sequence almost 
certainly contains a boarding group whose configuration 
is exactly front-to-back, e.g., {1, 2, . . . , m — 1, to}. Since 
boarding time for this particular group cannot be less 
than TO, so is boarding time of the entire passenger se- 
quence, i.e., limjv-^ooT' > to. With a lower bound of T 
thus obtained, we also build an argument for its upper 
bound. Let us choose one passenger from the sequence. 
When the boarding begins, the passenger has at most 
m + n< — 1 sites to go before reaching its seat, where 
n< denotes the number of passengers in the same group 
whose seat indices are smaller. Since those n< passengers 
are closer to their seats, they would get seated before the 
chosen passenger, leaving empty sites in the aisle. The 
chosen passenger advances along the aisle by at least one 
site per unit time, and can move further if empty sites ap- 
pear on its way. Hence, the chosen passenger can reach 
the reserved seat within to — 1 time steps and spends 
one more time step to sit down. Since the choice of the 
passenger was arbitrary, individual boarding time for any 



passenger is not greater than to, and so is the total board- 
ing time, i.e., T <m. Combining both upper and lower 
bounds, we obtain lim^r-foo T — m. 

Before moving on to the multi-column case, we remark 
that dividing into smaller groups always reduces board- 
ing time if c = 1, as previously reported by [16]. This 
is the case when keeping the passenger sequence ordered 
is always good for saving time. However, we emphasize 
that the benefits of group division, in terms of scaling, are 
not very robust against disorder caused by arrival time 
fluctuations. A slightest hint of globalized disorder can 
revert the scaling behavior to that of the random board- 
ing policy. Hence, division into small boarding groups 
is not very beneflcial, if one aims to reduce the scaling 
exponent of boarding time. 

Multi-column airplane — As a next step, we consider 
a fully reserved plane with multiple columns (c > 1) of 
seats, which is a very typical situation encountered in 
reality. Without loss of generality, we only focus on the 
boarding process of a two-column (c = 2) plane since the 
scaling behaviors does not change in the other cases. 

Figure [3]Je)-(h) show the scaling behaviors of (T) at 
c = 2. If no shuflling is applied, (T) scales as N^/'^ (N) if 
N/m (to) is fixed (see Fig.[He)). Gaussian shuffling at fl- 
nite TO results in exactly the same scalings (see Fig. [Sjf ) 
and (g)), implying (T) ~ TV^/^ ((T) - N) for purely 
globalized (localized) disorder. If both kinds of disorder 
are mixed together, the effect of globalized disorder pre- 
dominates, as the linear scaling for p < 1 indicates (see 
Fig.IHh)). 

For the special case when no shuffling is involved, the 
origin of the linear scaling can be explained as follows. 
Since the probability of a boarding group having the 
configuration like {1, 2, . . . , m, 1, 2, . . . , to} is finite, the 
number of such boarding groups grows linearly with N. 
When the first half of such a group begin to take seats, 
its later half completely occupy the rows reserved for the 
next group. This means that the next group cannot reach 
their seats for at least one time step, which implies that T 
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FIG. 4: (Color online) Asymptotic scaling behaviors of (T) 
in a partially reserved two-column (c — 2) airplane with 
group size m — 5. For the curves from top to bottom, (a) 
(p = 0.25,0.45,0.5,0.55,0.75 without shufHing (inset: semi- 
log plots of the lowest two curves in the main) and (b) 
(/) = 0.25, 0.35, 0.5, 0.55 with uniform shuffling p = 0.5. Each 
data point is averaged over at least 100 samples. 



increases by one or more time steps for every single group 
having this particular configuration. Thus, T > aN as 
N ^ oo, where a is a constant given by the probability 
of the configuration. Since T < N is always true, the 
leading order term of (T) must be linear in N. 

We note that division into smaller groups always in- 
creases boarding time if c > 1. In this case, reducing the 
range of disorder always results in stronger interference 
between passengers to be seated near each other, as re- 
flected in the strong sensitivity of the scaling exponent 
to localized disorder. Thus, lack of control is better than 
row-wise division of groups for c > 1 planes (provided 
that we do not divide passengers column-wise, which was 
covered by the c = 1 case), which is a lesson shared by 
most of previous studies [5l4l4 1 . 

Vacancy effect — At last, we consider the effect of va- 
cant seats on asymptotic scaling behaviors of (T) for a 
c = 2 case (see Fig. HJ, where each seat can be vacant 
with probability The same effect was also studied by 
01 , but only at a fixed plane size. 

Figure SJ^a) shows the effect of when there are no 
arrival time fluctuations. (T) scales as iV for (/) < 1/2 
and as IniV for > 1/2, implying that (f) can interpolate 
between the scalings for localized disorder observed at 
c = 1 and c > 1. The same observation can be made 
even in the presence of globalized disorder, as shown in 
Fig. Iljb), where the scaling changes from (T) ~ iV to 
(T) ^ N^/^. The transition point of the scaling exponent 
is 1/2 for c = 2, which is confirmed to be 0c = 1 ^ 1/c 
for the general value of c (even including c — 1, in which 
case the scaling is not affected by vacancy) . Interestingly, 
(T) ^ N^/^ seems to hold exactly at = </)c, regardless 
of the nature of sequential disorder. 

All the scaling behaviors of boarding time can be sum- 
marized as in Table m In a single-column (c = 1) plane, 
purely localized disorder produces the N'^ scaling, which 
includes both saturation and logarithmic divergence, and 



TABLE I: Summary of asymptotic scalings of (T). Note that 
Af" stands for both saturation to a finite value and logarithmic 
scaling. The gray arrow indicates that vacancy effect changes 
the multi-column scaling to the single-column one at ^c. 



Type of disorder 


Number of columns 
Single (c = 1) Multiple (c > 1) 


Localized 

(p = 0, cr ~ N^) 

Mixed 
(0 <p < 1) 

Globalized 
{p= 1,0-^ ATI) 


TVO/ 1 

/ Vacancy (f>: 

( Transition at 

\ <^c = l-l/c 
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purely globalized disorder leads to the N^^'^ scaling. The 
boarding process is more sensitive to globalized disorder, 
so the iV^/^ scaling is observed when both kinds of dis- 
order are present. On the other hand, in a multi-column 
(c > 1) plane, purely localized disorder produces the lin- 
ear scaling, while purely globalized one again yields the 
iV^/^ scaling. Since localized disorder is dominant in this 
case, the linear scaling is observed when both kinds of 
disorder are present. Finally, increasing the probabil- 
ity of vacant seats cj) beyond (jjc = 1 — 1/c changes the 
scaling behaviors from the multi-column ones to the cor- 
responding single-column ones while keeping the nature 
of sequential disorder. However, the borderline scaling at 
(j)c seems to be N^^^, regardless of the nature of disorder. 

To sum up, we have systematically investigated the re- 
lationship between the range of disorder in the passenger 
sequence and the asymptotic scaling behavior of boarding 
time using a simple boarding model. Our results clarify 
the origins of different scalings and indicate which type 
of disorder plays a dominant role. This offers a natural 
way to incorporate the boarding problem, previously re- 
garded as an engineer's optimization problem, into the 
domain of physics. 

Finally, we add a few remarks on how the results of our 
study can be further generalized. For example, we can al- 
low fluctuations in the time required for each passenger to 
sit down. Since the scaling behavior of the longest block- 
ing subsequence does not change, such fluctuations would 
not affect the scaling behaviors if they are finite and un- 
correlated, but diverging or correlated fiuctuations may 
produce interesting changes. Alternatively, we can gen- 
eralize the range of sequential disorder by applying Gaus- 
sian shuffling with a ~ A^* for < s < 1, which seems to 
produce intermediate values of exponents [2^. We also 
note that the model considered in our study is similar to 
asymmetric simple exclusion process (ASEP) [^H, which 
naturally raises the question of how the concepts of ASEP 
can be applied to the case of boarding. We hope these 
points to be satisfactorily addressed in future studies. 
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